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We review the determinations of the pseudoscalar glueball and eventual radial excitation of the r] 
masses and decay constants from QCD spectral sum rules (QSSR). The glueball mass is (2.05 ± 
0.19) which one can compare with the eventual experimental candidate X(1835), while the 77(1400) 
is likely a radial excitation of the r;'-meson. Their effects on the estimates of U{V)a topological 
susceptibility and its slope as well as the impact of the latter in the estimate of the spin of the proton 
is discussed. We predict the singlet polarized parton distributions to be a^{Q^ — 4 GeV^)=0.32 ± 
0.02, which is about a factor two smaller than the OZl value, but comparable with the COMPASS 
measurement of 0.24 ± 0.02. 



Prologue 

It is a great honour and pleasure for me to write this review for Professor Adriano Di Gi- 
acomo for celebrating his 70th birthday. My scientific connection with Adriano started 
when I studied at CERN in 82, the SVZ-expansion of QSSR^'^, and where the Pisa group 
^ has presented the first quenched lattice results values of the gluon condensates (a^G^) 
and fabc{G"'G^G'^) at the 1st Montpellier Conference on Non-pertubative methods in 
85. About the same time, my previous interest^ on the topological susceptibility x(0) of 
the U{1)a anomaly^ and its slope continued after several discussions with Gabriele 
Veneziano during my stay at CERN . A program on the estimate of these quantities using 
QSSR has started ^. In parallel, the Pisa group has runned their lattice simulations, and, 
for his 2nd participation at the Montpellier conference (QCD 90), Adriano has presented 
the lattice Pisa group^ results for x(0) and x'(0) which confirmed our previous results in 
Since then, there was a common interest between the Pisa group and Montpellier on these 
non-perturbative aspects of QCD. More recently, these common interests concern the proton 
spin problem^'^'^'^^, the estimate of the gluon condensate using 

QSSR 12 and lattice", and 
the estimate of the quark-gluon mixed condensates using QSSR^'* and field correlators . 
In 94, Adriano joined the international organizing committee of the QCD-Montpellier Se- 
ries of Conferences and, in 2001, the one of the Madagascar High-Energy Physics (HEP- 
MAD) Series of Conferences. He is among the few committee members who send regularly 
speakers and participate continuously to these conferences. The organization team of these 
meetings has always appreciated his human kindness and friendship. 



1 Introduction 

The U anomaly is one the most fundamental and fascinating problem of QCD, which 
has been solved in the large Nc limit by 't Hooft-Witten- Veneziano where Veneziano 

"This review is written in honour of Professsor Adriano Di Giacomo for his 70th birthday 
(Adrianofest 26-27th January 2006, Pisa, Italy) 



found a solution without the uses of instantons. The U{1)a topological susceptibiUty is 
defined as: 



where: 



X(0) = ^gg{0) = i J d'x{TQ{x)QH0)) , (1) 



Q{x) = J-a.Grc;^, , (2) 



Stt 

is the f/(l)A anomaly of the singlet axial-current: 



d^'Ai.ix) = J2 ['^1 = ^^qMi^5)'^q] + "^rifQix) . (3) 

q=u,d,s 

where n/ is the number of light quark flavours, ipq the quark fields, niq the current quark 
mass and: G'^^, = {ll2)e^vpaG'^ . The topological charge is defined as: 

Q = J d'^x Q{x) , (4) 

which is an integer for classical field configurations (2nd Chem number). Arguments based 
on large Nc for SU{Nc) x SU{nf) lead to^ : 

X(0) = ^ + m2 - 2Mi) . (5) 

where, = (92.42 ± 0.26) MeV* . Forn/ = 3: 

X(0) ~ (180 MeV)^ . (6) 

Gabriele Veneziano and, later on, Giancarlo Rossi have challenged us to analyze the validity 

of this approximation at finite Nc = 3, by computing x(0) and x'(0) using QSSR, because 
the previous result is obtained at = M^, but not at = 0, where, they have used the 
expansion: 

X(g') x(0) + q''x'{0) + (7) 



2 The pseudoscalar gluonium/gluebaU mass from QSSR 

QSSR has been used to calculate the topological charge in pure Yang-Mills QCD. In so do- 
ing, one has first extracted the mass and decay constant of the pseudoscalar glueball/gluonium 
using either a numerical fitting procedure^'^ or sum rule optimization criteria^^'^^ from the 
unsubtracted sum rules (USR): 

£{t) = f^' dte-^^lm xit) and 7^(T) = --^ log/:(r) , (8) 
Jo dr 

''A generalization of the result including SU (3) breakings for the decay constants can be found in^^. 



where the QCD model coming from the discontinuity of the QCD graphs has been used for 
the continuum. The analysis gives: 



Mg = (2.05 ± 0.19) GeV and = (8 - 17) MeV , (9) 
corresponding to = (6 — 7) GeV^, and where fo is normaUzed as: 

{0\Q{x)\G) = V2fGM^ , (10) 

i.e. like U/{2nf) with = (92.42 ± 0.26)(92.42 ± 0.26) MeV. The positivity of the 
spectral function gives: 

Mg < (2.34 ± 0.42) GeV . (11) 
One can compare the previous values with the quenched lattice results ^^: 

Mg = (2.1 - 2.5) GeV . (12) 

Recent sum rule analysis using instanton Uquid model leads to similar results: Mg — 2.2 
GeV and /g — 17 MeV in our normalization^*^, though the same approach leads to some 
inconsistencies in the scalar channel One should note that in the previous analysis ^^'^^ 
and^*^, the value of the gluonium decay constant is smaller than the corresponding value 
of the ri' one which is about 24 MeV in the chiral hmit^'^'^ (see section 5.2). In a recent 
paper ^ , the mass of the eventual experimental pseudoscalar gluonium candidate X{1835) 
has been used as input and the Laplace sum rule C has been exploited for fixing the corre- 
sponding decay constant. Though (a priori) self-consistent, this analysis is less constrained 
than the previous sum rules used in^^'^°, where the two sum rules C and TZ have been si- 
multaneously used for constraining the gluonium mass and decay constant. The resulting 
value of the gluonium decay constant is about (8-12) MeV^^ and agrees with the previ- 
ous values. On the other, a G-rjc mixing due to direct instanton has been considered in 
for pushing the unmixed gluonium mass of about 2.1 GeV down to 1.8 GeV with a mix- 
ing angle of about 17°. This value is comparable with the early value of about 9p = 12° 
of the meson-gluonium mixing angle using the OPE of the off-diagonal light quark-gluon 
correlator: 

^ggiq') ^ 2^ / d^xe^'^^TQixplm , (13) 

which is proportional to m^^^. Indeed, it is also plausible that the X(1845) comes from a 
comphcated mixing of the glueball with the 77(1440) and ijc- A confirmation of the gluonium 
nature of the X(1835) requires some more independent tests. 

3 On the natures of the r?(1295) and 7?(1400) from QSSR 

There are two other experimental candidates which are the r?(1295) and ry(1400) which 
one can intuitively interpret as the first radial excitations of the ry(547) and ry'(958). Here, 
we shall test if the ?7(1400) satisfies this interpretation. Using the same approach, the nature 
of the 77(1 400) has been tested by measuring its coupling to the U(1)a singlet current 
Q{x) in the chiral Umit. In this case, we can work with the SSR: 



f*- dt 1 ^ , , 
/ — e *^-Imx(t) , 

Jo t IT 



(14) 



as one expects from large Nc arguments that the topological charge x(0) vanishes due to the 
6'-independence of the QCD Lagrangian, thanks to the presence of the singlet rji. Including 
in the sum rule, the contributions of the r]' , 77(1400), G and the QCD continuum with tc 
fixed previously, one cannot find any room to put the 7/(1400), i.e. k, 0. Relaxing the 
constraint by replacing the QCD continuum with the 7/(1400), one can deduce: 

/r,(i4oo) < 16 MeV < ~ (24.1 ± 3.5) MeV , (15) 

which should be a weak upper bound within this assumption. This feature indicates that 
the 7/(1400) is Ukely the first radial excitation of the t/', as intuitively expected, while the 
glueball has a higher mass. 

4 The U{1)a topological susceptibility from QSSR 

4. 1 The pure Yang-Mills result 

Applying the Borel/Laplace operator to the subtracted quantity: 

X{Q ) ~ X(Q)|no quarks i\6) 

one can derive a combination of the unsubtracted sum rule (USR) and subtracted sum rule 
(SSR) for the topological suceptibihty in pure Yang-Mills without quarks ^'^: 

X(0)|no quarks = fj J C''^ [l - ^-^^ hv^ x{t) - 

©'l^-lsiW + ^-^^^^'+O-'^^^V}. (.7) 

The value of the gluon {agG"^) = (0.07 it 0.01) GeV^ from e+e^ data and heavy-quark 
mass-splitting ^^'^ has been confirmed by lattice calculations with dynamical fermions 
while the value of the triple gluon condensate fabdG'^G^G'^) is about 1.5 GQV'^{asG'^) 
from instanton model and lattice calculations ^ . The previous combination of sum rules 
is more interesting than the alone SSR, as the effect of the QCD continuum is minimized 
here. At the sum rule optimization scale of r about 0.5 GeV""^, it leads to the value ^'^: 

X(0)|no quacks - -[(106 - 122) MeV]^ « -4 (^)' (G^) , (18) 

indicating the role of the gluon condensate in the determination of x(0). The sign and the 
size of this result, though inaccurate are in agreement with the large Nc results obtained 
previously. In pure gauge, lattice gives the value ^'2^'^^: 

xmsu{2) ^ -[(167 ± 25) MeMf , x(0)|5C/(3) - "[(191 ± 5) MeV]^ , (19) 



which one can compare with the two former ones from large Nc and from the sum rules. 



4.2 Result in the presence of quarks 



In this case the analysis is more involved as one also has to consider the diagonal quark- 
quark correlator'^: 

^ j d^x e'^^TJ,{x)Jlm , (20) 

and off-diagonal quark-gluon correlator in Eq. (13). Then, the full correlator reads: 

^5(9') = i^ggiq^) + 2V'9,(?') +%ci{q^) . (21) 
In this case, ^5 (0) is not vanishing and can be deduced from chiral Ward identities to be^'^°: 

V'5(0) = --7^^ "^^(w) , (22) 

^ q=u,d,s 

Therefore, the topological charge x(0) vanishes in the chiral limit due to the 6 independence 
of the QCD Lagrangian. Lattice calculations in full QCD for two degenerate dynamical 
fermions find^^: 

V'sW = [(163±6)MeV]S (23) 

in agreement with the large Nc result ^'^ but does not have enough accuracy for checking 

the linear mg-dependence expected from current algebra. 

In this case of massive quarks, we include SU (3) breakings to the sum rule in order to see 
the effect of the physical rj'. Including both the gluonium and QCD continuum contribution 
to the appropriate sum rule, one can see from'''^ that their contributions tend to compensate, 
and then give the approximate numerical formula to leading order: 

.MJ/K«>(:-^).,.(„).J,(I|)%-. (24) 

Using the quark model prediction: 

Z^, ~-^V3/^~27MeV, (25) 
2n/ 

the physical r]' mass and the value rns{r) ~ 100 MeV"^", we can deduce at the sum rule 
optimization scale r 2± 0.5 GeV~^: 

V'5(0) ~ [157MeV]^ , (26) 

in good agreement with the large Nc and lattice results, then showing the consistency of 
the sum rule approach. Armed with these consistency checks, we shall now evaluate with 
confidence the slope of the topological susceptibiUty x'(0). 

"These correlators include in their definitions the quark mass through Jq, while in^°, this quark mass is 
factorized out. 



5 The slope x'(0) of the topological susceptibUity from QSSR 

Using the expansion in Eq. (7), one can derive a sum rule for x'(0) by applying, the 
Laplace sum rule operator, to the twice subtracted quantity: 



x(g')-x(o)-gV(o) 



2^2 



{1 

5.1 The pure Yang-Mills case 
One obtains in this case: 

dt 



[^-^--Iraxit], (27) 

J t t — — tSTT 



f 

Jo 



*^-Imx(t) = ^i(r) - x(0)t + x'(0) , (28) 



where ^i(r) comes from the OPE expression of the correlator x(?^)^- Eliminating x(0) 
by using Eq. (17), one can deduce the sum rule^: 



oo 



dt 



e *^ 




l-tr(l-l) ilmx(t)]:.x'(0)+fe)'^r-x 



|l+ (^^) [lO-27BA-2^1og(-logrA2)] + ATr'^giG'^^ , (29) 

where Pi = —11/2 and P2 = —51/4 for pure gauge SU (3)c. Using r- and ic-stabiUty cri- 
teria, and saturating the spectral function by the pseudoscalar glueball and QCD continuum, 
one can deduce at r ~ 0.5 GeV~^^: 

X'(0)|no quarks - -[7 ± 3) MeV]^ , (30) 

indicating that: x'(0)M^, <^ x(0) ) which justifies the accuracy of the large Nc result. This 
value and the sign has been confirmed by lattice calculations^: 

X'(0)|quenched ^ -[(9-8 ± 0.9) MeV]' . (31) 
5.2 Result in the presence of massless quarks 

In this case, quark loops enter into the value of the j3 function and of the QCD scale A. The 
■q' enters now into the spectral function. Its coupling to the gluonic current can be estimated 
from the SSR: 

^e-^^fl-^l-ImxW, (32) 



taking into account the fact that in the chiral limit x(0) = 0. Therefore, one obtains^: 

/^/ (r ~ 0.5 GeV-2) ~ (24.1 ± 3.5) MeV , (33) 
where the decay constant has a weak r-dependence because of renormalization: 

/,,(r)^/,,exp(-^^A__) , (34) 



where f^^ is RG invariant. This result can be compared with the quark model prediction 
in Eq. (25). Using the Laplace and FESR versions of the twice subtracted sum rule, one 
finds^: 

Vx^lma.siess(r) 0^ [(26.5 ± 3.1) MeV]' . (35) 

The result has changed sign compared to the one from pure Yang-Mills, while its absolute 
value is about 12 times higher. The main effect is due to the r]' which gives the dominant 
contribution to the spectral function compared to the gluonium and QCD continuum. That 
can be understood because of its low mass and of the fact that its decay constant is larger 
than the gluonium decay constant fc- 

Similar estimate of V'5 (0) in the chiral limit for the singlet channel has been done in^^ , where 
the value of x(0)' is significantly larger by about a factor 2.5 than ours. However, we found 
some inconsistencies in the two sides of the SR: in the experimental side the contributions 
of the singlet and octet mesons have been considered and the phenomenological value of 
the decays constants, masses and mixing angles have been used; while in the QCD side, 
only the correlator associated to the singlet current Q{x), or to massless quark has been 
accounted for. 

On the other, some criticisms raised in a series of paper due to direct instanton breaking of 
the OPE, do not apply in our case: our different sum rules optimize at a large scale r < 0.5 
GeV~^, where this contribution become irrelevant like other high-dimension condensates. 
In^'', it has been argued that screening corrections cancel the direct instanton contributions, 
indicating that these effects are not yet well understood. However, in your approach, the 
introduction of the new l/tj'^-term in the OPE, which is also negligible, is an alternative to 
the direct instanton contributions Some answers to these criticisms have been akeady 
explicitly given in^*^. 



5.3 Result in the presence of massive quarks 

As can be seen from the expression of the two-point correlator VsCq^) in Eq. (21), the 
analysis is more involved. In this case, its slope has been estimated from the substracted 
Laplace sum rules. At the stability point r c± (0.2 — 0.4) GeV~^ , one finds 

y^[(0) = (33.5 ± 3.9) MeV , (36) 

while the r/' -decay constant becomes: 

f^, (r) ~ (27.4 ± 3.7) MeV , (37) 

compared with the results for massless quarks in Eqs. (35) and (33) and agrees quite well 
with the quark model prediction in Eq. (25). One can notice that the SU (3) breaking effect 
is about 10% and 20% respectively showing a smooth dependence on the strange quark 
mass as expected. Similar analysis can be done in the flavour non-singlet case by working 
with the correlator associated to the 77-meson current: 

(0|9^J^5|r/) = /,M^^ (38) 

where in the SU{3>) hmit and in this normalization without l/2nj, /»j = /tt = (92.42 it 
0.26) MeV. Therefore, one obtains ^°: 



^Ji^^iG) = (43.8 ± 5.0) MeV , and /^//^ = 1.37 ± 0.16 



(39) 



Instead of the value of /^^p = (92.42 ± 0.26) MeV, we have used the sum rule prediction 



10. 



U = (107 ± 12) MeV , (40) 

for a self-consistency of the whole results. The value of the SU{3) breaking ratio fri/ fir is 
in Une with fx I fix = 1-2, where we expect bigger effects for the rj than for the K. 

6 Applications of the QSSR results to the proton spin 

The previous results have been applied to the proton spin problem, where one expects that 
the gluon content of the proton is due to the ^(1)^ anomaly. This property can be made 
explicit in the approach of DIS where the matrix elements from the OPE are factorised into 
composite operators and proper vertex functions . In the case of polarised fip scattering, 
the composite operator can be identified with the slope x'(0) of the topological suscepti- 
biUty, which is an universal quantity and then target independent, while the corresponding 
proper vertex is renormalisation group invariant. The first moment of the polarised structure 
function reads, in terms of the axial charges of the proton: 

(41) 

where the Wilson coefficients arise from the OPE of the two electromagnetic currents: 

C^^ = l-as- 3.583a^ - 20.215af , Cf = 1 - Ja^ - 0.5500^ , (42) 

where Ug = as/n. The axial charge are defined from the forward matrix elements as: 

{P, s\J^5\P^ s) = ^a^s^ , {p, s| J^sIp, s) = ^^(^^s^ , {p, s| J^sIp, s) = aP{Q^)s^ , 

(43) 

where J^g are the axial currents and the proton polarisation vector. Using QCD parton 
model, the axial charges read, in terms of moments of parton distributions^^: 

a^ = Au-Ad, = Au + Ad-2As , a° = [AE = Au + Ad+As] -^OsAc/(Q^) . 

(44) 

and are known in terms of the F and D coefficients from beta and hyperon decays: 

= F + D and = 3F - D , (45) 

where ^^: 

F + D = 1.257 ± 0.008 and F/D = 0.575 ± 0.016 , (46) 

so that an experimental determination of the first moment of in polarised deep inelastic 
scattering (DIS) allows a determination of the singlet axial charge aP{Q'^). In the naive 
quark or valence quark, parton model, one expects that As = Ag = 0, and then a° = a®, 
which is the OZI prediction. The proton spin problem is that the experimental value of 



is much smaller than a^, which would be its expected value if the OZI rule were 
exact in this channel. The first SMC data found ^^'^^ : 

a°(Q2 = 10 GeV^) =0.19 ±0.17 and o"(Q2 = 5 GeV^) = 0.19 ± 0.06 , (47) 

which are confirmed and improved by the recent COMPASS data at = 4 GeV^^^: 

aO(Q2 = 4GcV^) = 0.24±0.02 . (48) 

These results are much smaller than the OZI prediction {Ag = 0): 

AS|ozi = 3F-Z) = 0.579 ±0.021. (49) 

Inserted into Eq. (41), the OZI results lead to the Elhs-Jaffe sum rule It has been 
conjectured in^^ that the suppression of a'^ with respect to is due to the gluon distribution. 
However, the SMC measurement of AE gives a value 

AS|sMC- 0.38 ±0.04, (50) 

indicating that independently of Ag, there is also a large difference between the OZI predic- 
tion and the data'^. In the following, we will not try to solve this discrepancy but will show 
our prediction for a°, where the sum of the quark and gluon components are concerned. 

6.1 Prediction in the chiral limit 

Using the composite operator © proper vertex functions approach^, one can write the de- 
composition, in the chiral limit: 

r? singlet = l^'2rifCf{as{Q'))./^)\Q2t^or,M . (51) 

The key assumption is that the vertex is well approximated by its OZI value: 



while all OZI violation in singlet contained inside ^/x^Wj. Comparing the result with 
the OZI prediction of a^, one can deduce: 

= ^\f^)\Q^ = 0-60 ± 0.12 , (53) 

which gives our original proton spin sum rule^: 

a°{Q^ = 10 GeV^) = 0.35 ± 0.05 =^ r?(Q2 = 10 GeV^) = 0.143 ± 0.005 , (54) 
where we have used: 



=10 Gev2 = (23.2 ± 2.4) MeV , (55) 

after running the result in Eq. (35) from 2 to 10 GeV^. These results can be compared with 
the OZI value in Eq. (49) and agree with the last SMC data at = 10 GeV^^°: 

rf(Q2) = 0.145 ± 0.008 ± 0.011 =^ a°(Q2) = 0.37 ± 0.07 ± 0.10 , (56) 



'^In*^, an estimate of the scalar sea quark content of the nucleon indicates that it is suppressed compared to 
the valence one. Similar results may apply here. 



6.2 The case of massive quarks 

In this case, the relation in Eq. (53) is replaced by^°: 

' = ] = 0.55 ± 0.02 , (57) 

" V2y>r(0) 

where ip5{q'^) has been defined in Eq. (21) and tp^^iq^) is the two-point correlator for the 
octet current. The running of the subtraction constant ipi,{0) is very smooth from r~^= 3 to 
10 GeV^. Using V/^^(0) = = 0.58, the previous relation leads to: 

o°(Q2 = 4 GeV^) ~ a°(Q2 = 10 GeV^) ~ 0.32 ± 0.02 , (58) 

which is almost equal to its value in the chiral limit. This prediction is comparable with the 
COMPASS result (0.24 ± 0.02)^^ given in Eq. (48) but is about a factor 2 smaller than its 
OZI value in Eq. (49). Using this result, we predict: 

(Q2 = 4 GeV^) = 0.145 ± 0.002 , (59) 

where we have used the value of as(M,-) ~ 0.347 =b 0.03^^'^"^'^^, and the previous values 
of F and D from Improvements of these results at COMPASS energy require a good 
control of the higher twist corrections in the PT expressions of which are expected to be 
bigger here than at SMC. Using the sum rule estimate in^^ of about —0.03 of the coefficient 
of the 1st power corrections, one obtains a correction of about -(0.008 it 0.008) to 
at 4 GeV2, where the error is an educated guess taking into account larger absolute 
value obtained from the fit of the Bjorken sum rule in^. A control of the corrections in the 
assumption of the validity of OZI for the singlet and non-singlet vertex functions is not yet 
testable^. Including the power correction, we consider as a final estimate: 

rP{Q^ = 4 GeV^) = 0.137 ± 0.008 . (60) 

Several tests of the approach have been also proposed in the Uterature"^''. Another crucial 
test of our result will be a lattice measurement of (0) with massive dynamical fermions, 
which we wish that the Pisa group puts in its agenda. Unfortunately, an attempt to measure 
aP{Q'^) on the lattice has not been conclusive, as it gives a value: 

aO(Q2) = 0.04 ±0.04 ±0.20, (61) 
where the last error is an educated guess of the effective error expected by the authors. 



7 Conclusions 

We have used QSSR for predicting the gluonium decay constant and mass, and for arguing 
that the 77(1440) is likely the radial excitation of the ?7'(958). These results have been used 
for predicting the topological susceptibility and its slope, which are useful inputs in the 
resolution of the "proton spin problem" . 



"A QSSR evaluation of the ijNN coupling has not been conclusive 
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